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^ Abstract 

^~~^ A model of the charge-transfer process in DNA is investigated by using the Feynman path integral 

^^ [m [5T1 [321 [Sn] . The base-pairs' coordinates are eliminated resulting to the effective action of the system. 

. The trial action is introduced in order to obtain the expectation value of the equation of motion. The 

^^ electron ■will be in the equilibrium bet^ween lose and gain energy under the presence of an external electric 




field -when the steady-state is reached. Using the result of the steady-state condition, the impedance 
function of the charge moving along the DNA chain is also studied. The mobility of charge moving in 
• • DNA ■was obtained and the limits at lo'w and high temperature are also studied. The results of ho'w 

. 5^ the mobility depends on temperatures indicate that DNA may possibly be treated as a semiconductor 

S^ qunatum wire. 

J3 



1 Introduction 

DNA is one of the most important and complex structures of living organisms. In a nut shell, DNA 
is a nucleic acid containing the genetic instructions and plays a major role of the long-term storage of 
information, or a set of blueprints. The structure of DNA had been found to be the double helix, by 
Watson and Crick [5S] at Cambridge University in 1953, composing of two polymers. These two polymers 
are glued together by base pairs which are the Adenin(A)-Thymine(T) and the Guanine(G)-Cytosine(C) 
pairs, see Fig. [T] Understanding all possible physical properties of DNA is one of the main quests of 
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Figure 1: DNA is a double helix connected by the base-pairs, i.e., A-T and G-C [65] . 

Biophysicists. Studying the nature of charge transport mechanism in DNA molecule may gain more basic 
knowledge about signaling and repairing process of damaged DNA that are related to the mutation of 
the DNA resulting to cancer. The mechanism of charge transport through DNA could also be able to 
use to protect the damage, which is caused by oxidizing chemical [45l [20] , of genomes in some organisms. 
Moreover, DNA could be possibly used as a tool to develop numerous tiny devices such as DNA biosensors, 
template for assembling nanocircuits and electric DNA sequencings, which are expected to be different 
from those existing traditional electronic devices [121 [331 \ST[ [551 [33] . 

There were many experimental results which indicate that DNA could be a conductor. The first 
confirmation was from Eley and Spivey ^ initiating that DNA exhibits conductor property. Consequently, 
scientists have been attracted to investigate the electronic properties of DNA. DNA could be a conductor 
because of the formation of a 7r-bond across the different A-T and C-G base-pairs. The 7r-stacking 
interaction is a key function for the electrical properties of other organic molecules. Fink et al., [Ill 
[T2] have shown the result of I-V characteristic of DNA suggesting that DNA is a conductor. Later, 
Porath et al., p3] have also been created an experiment to show that DNA could be treated as a large- 
band gap semiconductor. However, the charge transfer mechanisms in DNA chain are not yet well 
understood and remain very controversiaF] At temperature above 1 K, DNA may behave as conductors, 
semiconductors, insulator, while at below this temperature, DNA may exhibit a property as proximity- 
induced superconductor [321 [Ml El [331 [SI] ■ Interestingly in ^T] , the electric conductivity of A-DNA 
using DNA covalcntly bonded to Au electrodes was investigated. They found that A-DNA is an excellent 
insulator on the micron scale at room temperature. 

Many theoretical models proposing to study the electrical properties of DNA were based on incoherent 
phonon-assisted hopping [321 I25j , classical diffusion under the conditions of temperature-driven fluctua- 
tions i31, variable range hopping between localised states 57|, transport via coherent tunneling [S' and 
polaron and soliton mechanisms jSi i33t i22i i23t ii40i, i56j . The quantum model of charge mobility has been 
studied by Lakhno (SO] [31] . The Kobu theory was needed to approximated the model yielding that the 
mobility increases with the temperature decrease. Another approach is Feynman path integrals qJ which 
was used to study the ground-state energy and the effective mass of the electron in DNA by Natda 39 . In 
the present paper, we will concentrate on the theoretical description of the mobility of the charge moving 

^The literatures on DNA conductivity can be found in [lOi rfllQ] 

^The Feynman path integrals was successfully used to study the dynamic of the electron in the polar crystal |13lil4|[T^I5UI 
[4Tl[21][4i[58l|59l[io] 



along a model of one- dimensional DNA chain. We use the path integral as a key technique to study the 
problem. The Hamiltonian of the system can be represented in the analogue form with Frohlich's Haniil- 
tonian suggesting that Feynman approach in polaron problem is applicable. We will show both analytical 
and numerical results of the study and we will try to give evidences: Could DNA be a conductor and what 
are the conditions and factors effecting to its conductivity. The structure of the paper is organised in the 
following. In section [2] we will discuss the model of one-dimensional DNA. The path integral method is 
used to eliminated the coordinates of the harmonic chains in DNA leading to the effective action of the 
system. In section [31 the steady-state condition at finite temperatures of the charge in DNA chain under 
the presence of the external electric field is studied. In section [41 the impedance function is computed 
from the equation of motion by relaxing the steady-state condition. Furthermore, the mobility at low 
temperatures of the charge is studied with varying the coupling constants and number of base-pairs. The 
effective mass at zero temperature of the charge interacting with DNA chain are also calculated. In sec- 
tion [5] the physical concept of the analytical analysis and numerical results in the previous sections will 
be presented. In section [6j the summary of results and some open problems along with possible future 
developments will be mentioned. 

2 The model of charge transport in one dimensional DNA chain 

DNA composes of two major long polymers which are made from repeating units called nucleotides. For 
all living things, DNA is in the shape of a double helix which is bonded together with backbones made 
from alternating phosphate and sugar residues, i.e., A-T and C-G pairs. As mentioned that DNA has 
a double helix structure which is complicated form. Then in this paper, the structure of DNA will be 
simplified by stretching a DNA chain into one dimension as shown in Fig. [2] The aim of the current 
study is to show the effect of the length, coupling constant and temperatures to the charge mobility in 
DNA. The electron is now moving along the DNA chain with N base-pairs and there is the interaction 
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Figure 2: A model of one-dimensional DNA chain. The base-pairs, which will be treaded to have the same 
physical structure, are represented by oscillators which will move back and forth along the x-direction. 

between the electron and base-pairs throughout the chain. The Hamiltonian of the model proposed by 
Nat da [3^] is given by 
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(2.1) 



The first term is the kinetic energy where p and m are the momentum and mass of the charge, respectively. 
The second term represents the harmonic oscillators of the ith base-pair(backbone) where M and J7 are 
the mass and oscillation frequency for each base-pair and Xi is the displacement along the x-axis of the 
ith from its equilibrium xo,i. The last term describes the interactioirM which is modeled by Dirac delta 
function, between electron and the base-pairs where a is the coupling constant. 



We assume that the interaction will take place only on the ith base-pair(distorted), not for other base-pairs(undistorted). 



The action of the system is now 
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(2.2) 



and the propagator can be written in the form 
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Using Feynman path integral method [131 ESI 1151 [3S]i the base-pairs' coordinates can be ehminated 
resuhing to the effective propagator 
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(2.4) 



where the S^ff is the effective action given by 
Spff = dt 
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where fi{t, a) — S[x{t) — xo^i]S[x{a) — xo^i]. Using the fact that J2i S[x{t) — xo^i]S[x{a) — xq,,;] ~ NS[x{t) 
x{a)], the effective action becomes 
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since 5[x{t)~x{a)] — -^ /^ dfce*'^^^^*^ '^^'^^^ and we define 77 = 



MWa^N 



After ehminating the base-pairs' 



coordinate, the action of the system (2.6) can be written as a function of the electron coordinate. 



3 The steady-state condition 

In this section, we consider the system of the electron moving along the DNA chain under the influence 
of the external constant electric field Eq. The Hamiltonian of the system now is 
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When the electric field is turn on, the electron will be accelerated. However, electron will interact with 
the base-pairs while it moves along the DNA chain resulting to loose some amount of energy to base-pairs. 
Then if we leave the system long enough the system will reach the steady-state. 

In order to get the steady-state condition, we will apply the method, called Double path integral 
[131 Sni HOI [SD] J that was used to study the motion of the electron in the polar crystal. Let p be the 
density matrix of the system and the expectation value of any operator Q at time t given by 



<Q(t)>-Tr(Qp(t)) 



(3.2) 



We now assume that the system is in the thermal equilibrium initially, p(t = 0) = exp(— /SiJe-DAr^) where 
/? = l/fc^T, T and ks are the temperature and Boltzmann constant, respectively. To obtain the density 
matrix at any time t, we need to solve the time-evolution equation 



ifl-KT = [He-DNA,p] 



(3.3) 



and therefore we havq3 

p(t) = e-iJ'o^"^''>'''p{0)e+TiJ'o^"''^''^'^'' . (3.4) 

As we mentioned earlier that at the beginning the DNA is in thermal equilibrium, then we may write 
p(0) = 6xp(— /3^j, alak) where a|, and a^ are creation and annihilation operators of the base-pairs. 



Next, we will compute the expectation value of electron's displacement at time T 

id J 
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where we now put Eq -^ Eo+jS{t — T) and Eq —> Eq. The problem of integrating Trp(T) over DNA base- 
pairs' coordinate can be performed in the same fashion as we mentioned in section [21 For the convenient, 
we transform the electron coordinate to a frame of reference drifting [121 [TH \5U] with the expectation 
value of the electron's velocity v. x{t) — y{t) -\- vt. Then Trp(r) can be written in the form 

Trp(T) = J ^[y{t)] J ^[y'(t)]e'*= , (3.6) 

where 
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Here comes to the crucial point. The path integral of (3.6) could not be performed exactly, then we need 
to introduce a trail action given by 
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The function G(r2') is the distribution function 

particle as shown in Fig. [SJ We now expand (|3.6 1 in the form 



[l4l |49] with frequencies O' = V x^ of j th fictitious 



''An unprimed operator acts to the left and ordered right to left with increasing time, and primed operator acts to the right 
and ordered left to right with increasing time. 

^In reference [131 1141 1491 I5UI I39| . they considered only one fictitious particle. The reason why we consider A^' fictitious 
particles is to make trial action more realistic and also to improve the variational result [I] |2] . 




Figure 3: The electron is attached to N' fictitious particles of mass M' with different coupling constant Kj 
corresponding to each one of them. 
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Using (3.10), we possibly write 
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The velocity of the electron can also be expanded in the form 

■^e = ''^tr + e^tr ~ tr'^tr + 
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We are interested to calculate the expectation value of the velocity at T — > oo which the steady-state is 
reached. Then we have 
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From the relation (3.13), we may find 
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which is noting more than the first order expansion of the reciprocal of v^ — Vtr + eVtr— tr^tr ■ The reason 
why we choose such an expansion for this type of problem is given in reference [TU 1321 [SD] ■ Moreover, 
assuming vtr = — trVtr which reduces (3.15) to 
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From (3.12b), we find that 
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where 
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The first order of velocity in the expansion is 

vtr = ~i2EahnL{T) . (3.19) 

The next velocity in the expansion ^Vtr can be computed by considering the term J J 'fee**'' 
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We introduce the shorthand notations where the force terms J J[Eq, E'q] in J J ^^e**'' wih be inserted as 
follows 



(1) = J J [Eo+jS{t-T),E'o + k{S{t-e)-S{t-0)] , (3.21a) 

(2) = I l[E,+jS{t-T) + k{S{t-9)-S{t-0),E'o] , (3.21b) 

(3) = f f [Eo+-fS{t-T) + kS{t-e),E'„ + k5{t-0] , (3.21c) 

(4) = f f [Ea+j6{t-T) + k6{t-0,K + kSit-e)] . (3.21d) 



After performing a tedious computation, we put eft,, into (3.16) together with another assumption that 
Ve — Vtr = V yielding the steady-state equation 
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and Zj are the roots of the polynomial 
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(3.22) represents the equation of motion in the steady-state condition, which holds at any temperatures, 



for the charge moving along the DNA chain according to the model (2.1). The charge gains the energy 



from the electric field and gives the energy to the base-pairs via the interaction . Note that we compute 



the velocity only the first oder in the expansion (3.13) corresponding to the fact that the interaction 



between the electron and base-pairs is weak. The expansion in (3.11) can be pushed further in order to 
cover the higher range of the interaction. 

4 Electrical properties of the DNA 

In the previous section, the steady-state equation is obtained by looking the expansion of the velocity of 
the electron moving along DNA chain. For the present section, we will compute the impedance function 



of the system from the equation of motion. The electric field is now modified to Eq -> E{t) = E^e^^* 
where vj is the oscillation frequency. From (3.171, we may write 
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Taking the Fourier transform of (4.1 ), we have 

where Y^ is the Fourier transform of the response function given by 
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The Z^ is the impedance function of the electron moving on DNA chain. Without imposing the steady- 
state condition, the equation of motion reads 



/OO /•O 

dc 
-OO J — C 



smh[i2p/2J 



Applying (4.1 1 to (4.4 1, we find that the impedance function can be expressed in the form 
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We are now interested for the case of tn < 1 and j3 — )■ cx)(zero temperature) and the path of integration 
along the real axis may be rotated to a path along the positive or negative imaginary axis ^ ~ to ±ioo 
|14) . The impedance becomes 
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In the case of extremely low frequency nj, we may approximate 1 — coshfnj^] ^ ~vj^^ jl. Then the 
impedance function is simplified to 

Z^ = ~m^m* , (4.9) 



where m* is the effective mass given by 
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(4.9) together with (4.10) suggest that the charge behaves like a free particle with the effective mass m* 
when it moves along the DNA chain. According to the model, (4.7) is the impedance function for the 
whole DNA chain which consists of N base-pairs. If we assume that for each base-pair has the same 
impedance, we find 

^- = ^ , (4.11) 



which is the impedance for each base-pair. Then the system can be presented as the connected series of 
the impedance function shown in Fig|4j In the case that N becomes very large N ^ oo the impedance 
function remains finite 
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From the result (4.11), we see that the conductivity of the DNA is directly propor- 



where 77' 

tional to the reciprocal of the length of DNA j63j . 



Figure 4: DNA chain can be represented as the connected series of the impedance function z^ 



We now introduce the function 
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which is the last term of the impedance function ( |4.5[ ). For the convenient for analytical analysis as 
mentioned in [Ml |49l [50] , we may consider 
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The mobility [131 (Ml [SO] of the charge moving along DNA chain is defined 
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We see that the mobility of the charge is a function of the temperature and also the coupling constant a. 



We also would like to point that (4.141 and (4.161 hold at any temperatures. In the next section, we will 



show the analytical analysis as well as the numerical results of the effective mass and the charge mobility. 

5 Physical concept of the results 

In this section, we will focus more the physical meaning of the results from the previous sections. Before 



proceeding, we need to reduce the complexity of the trial action (3.8). Here, we are interested in the case 



of N' = 1. Let V and 51' be the variational variables, determined from the minimization of the ground 
state energy [33 [12] j given by 
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The steady-state condition: With N' = 1, the (3.22) becomes 
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We will follow the analysis process of (5.2) as those in [491 |4T]. We start to rewrite (5.3b I in the following 
form HI] 
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Inserting (5.5) and (5.4 1 into (5.2 1 together with expanding the terms [1 — n{v)\ and n{v) in the Taylor 
series, we obtain 
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The physical meaning of (^5.6| is the following, see also [491 HI]. The term in the second line of (5.6), with 



the term [1 + n{fl)], represents the process that the electron with momentum p + AIv emits a phonon 
with energy f2 and momentum k, while the electron state changes from n' — )■ n as shown in FigMa). 
The electron momentum is reduced by k and the delta function tells that energy is conserved, while the 



[1 + n{Q)] shows the a phonon with frequency Q is emitted. The term in the third line of (5.6), with the 



term n{fl), describes the opposite situation as follows: the electron with momentum p + Mv absorbs the 
phonon frequency 51 as shown in FigJMb). 




(a) 



p + Mv - k p + Mv 




Figure 5: The scattering processes of the electron (a solid line) by phonon (a wavy line), (a) a electron with momentum 
p + Mv in n' state emits a phonon with momentum k and energy fi. The electron changes to a new state n with 
momentum p + Mv — k. (b) shows the opposite situation in which the phonon is absorbed. 



The effective mass: The effective mass (4.10) can be simplified onto the form 
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Figure 6: Coupling interaction dependence of the effective mass for different number of base-pairs. The dashed, 
dotted, dot-dashed and two-dot-dashed lines are dependence of the effective mass when the number of base pairs are 
N = 30, 90, 120 and 1500, respectively. 
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which agrees with '39] (with replacement p with N) . The numerical results of the effective mass with vari- 
ous coupling constants and the number of base-pairs are shown in Fig. 6. In the case of the weak coupling, 
the electron will move freely as the effective mass becomes the mass of the electron: m* — > to = 1. In the 
case of the strong coupling, the electron will behave like a heavy particle as the effective mass increases. 
In the case of damaged DNA, the coupling constant of the base-pairs may possibly change to be in the 
strong coupling regime or another way around, then there is a potential to use the effective mass of the 
electron to detect the damage in DNA. Another effect to the mass of the electron is the number of the 
base-pairs N. For the large number N, the effective mass of the electron will increase drastically. This 
result suggests that the length of the DNA has effect to the charge transfer mechanism as pointed in [53] . 



The mobility: We consider the function Imxm for the case N' — 1 

cos[wu\ cos[u\du 



47277/33/2 / u x 3 sinh[n7/3/2l 
Imxro = I 



TT \n'J sinh[/3/2] 7o 



u^ + ^+P (^^fT ) coth[/?./2] ~ /3 (^4^) ,^^ 



[Pv/2] 



3/2 



(5.9) 



We see that the integrand of (5.9) has the same structure with those (47a) in [14,. If we now consider the 
system at low temperatures and ccj < 1 we have 



Imxc 



7r3^e-/3(i--)g- 



1(1-^) 



Vi 



zue ^e 






1(1+^) 



(5.10) 



The first term describes the process that the electron absorbs a quantum energy w from the electric field 
and emits a phonon with unit energy. The second term shows that the electron absorbs a quantum energy 
uj and also absorbs a phonon then at the end the electron will carry the momentum \/l -I- tjj. The third 
term, the electron absorbs a phonon and releases a quantum energy zu to the electric field. The fourth 
term represents the electron with energy \ + vj emitting a phonon as well as a quantum energy to the 

electric field. The term exp[— ( ^^72^ ) (1 ± ro)] tells the chance that the electron with energy 1 ± nj com- 
pletes the processes. Eq(5.10l contains possible scatterings happening while the electron moving along 
the DNA chain as shown in Fig. 7. 
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(c) 



(d) 



Figure 7: The scattering processes of tiie electron (a solid line) absorbing and emitting a quantum energy (a wavy 
line) and a phonon (a dashed line). 



In the case of low temperatures and tu -> 0, (5.9 1 leads to the dc mobility 

cos[u\du 



1 



4?7 



f^^"%ri: 



n I v^—iV^ \ 1 — cos[tJM] 



0'2 



P^ 



3/2 



(5.11) 



As we explained (5.101 that the electron may be scattered by oscillations of the base-pairs (phonons) 
which will be increased with the temperature. We are interested for the case of weak coupling because 
some furture analysis can be done on the mobility (5.11). In this case, the variational variables J7' « u 
leading to 

' 4,//?3/2 r cosMd. ^ ^^.llpy^Km; , (5.12) 



M 



.2 + f 



where Ki is the modified Bessel function of second kind. 

At high temperatures, /3 — > and Ki{f3) « 1//3, the mobility becomes 



1 1 

— oc —^ 



/i oc 



2.-1/2 _ 



(5.13) 



At high temperatures, electron-optical phonon interaction will play an important role in quantum wires 
of non-polar semiconductors which the mobility is affected strongly by the optical phonon scattering so 
the mobility is expected to be proportional to T~^/^ [B31[Tn]. For our result in (5.13), we observed that 
at high temperatures the behavior of the mobility in the case of DNA chain is similar to that in the semi- 
conductors quantum wires expression. We found n oc T~^''^, indicating that optical phonon scattering 
dominates the temperature dependence of the DNA mobility for weak coupling. 



At low temperatures, /3 — ?► oo and Ki{l3) « e^ /y/j3, the mobility becomes 



/z ex e 



/3 



(5.14) 



The electron mobility will be influenced strongly by the interaction of the electrons with phonons. In a 
semiconductor quantum wire, the mobility of an electron interacting with acoustic phonons would follow an 
exponential behavior and the energy relaxation of an electron due to emissions of acoustic phonon strongly 
increased by all low temperatures acoustic phonon scattering mechanisms, which are exponentially acti- 
vated in a an electron [5U [3j . In weak coupling systems of DNA chain, the low-temperature behavior of its 
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mobility has the characteristic form of an exponential term. The feature of our results is the identification 
of electron transfer in the DNA similar to be behavior of electron scattering in a semiconductor quantum 
wire. Recent studies of the electrical conductivity of DNA molecules reveal that they may act as semicon- 
ductor materials with nanometer scale dimensions [Sllig l [M l [55|[2^l?7 l l5 ^ I ^ I ^ I ^ ITT | [25 1 [^P71[7]. 

6 Summary 

We study the model of charge transfer in one-dimensional DNA chain with arbitrary iV base-pairs. The 
base-pairs have been modeled as the oscillators and the interaction between the charge and the oscillator 
is represented by the Dirac delta function. The path integral method is used as the key tool to study the 
physical properties of the system. The main advantage of path integral is that the base-pairs' coordinates 
can be integrated over leading to the effective action. The expectation value of the equation of motion 
of the electron moving in DNA chain under the influence of the external electric field can be computed 
by introducing the trial action. The electron will gain energy from the electric field and will give some 
to the phonons via the scattering process. In steady-state condition, the electron is the balance of these 
two processes. The electron will behave like a heavy particle if the coupling constant and the number of 
base-pairs are large. On the other hand, the electron will move freely in weak coupling region and small 
number of base-pairs. 

In the present paper, the mobility of charge moving along DNA chain is studied as a major feature. 
The study of mobility both at low and high temperatures suggests that DNA can be possibly treated as 
a semiconductor quantum wire. For the case of extremely low temperatures, the mobility is very high 
reflecting a very high conductivity. This result may be possibly explained by the resonant tunneling effect. 
According to our knowledge, the mobility of charge along DNA has not yet been measured directly and 
most of the experiments were focused at room temperature or high temperatures [9]. Then there is a 
difficulty to compare our results with the experiments, just some of them as we did explain in the main 
text before. 

The hopping process ^\ [TB] may be a current candidate to theoretically explain the charge transfer in 
DNA and indeed there are many experimental results supported |47ll55l[T8] . However, experiments have 
so many factors as well as conditions that need to be controlled leading to a wide range of the results. This 
causes a big debate on what should DNA conductivity be as pointed out in [7] . Then, we believe that 
the model in this paper may offer an alternative theoretical description to the charge transfer process. 

At the end, we would like to point that the model can be developed for more realistic. Then the 
results in the present paper can be improved by taking into account the following points. The first one 
is to consider the different coupling constants between A-T and C-G pairs. The second is to consider the 
anharmonic oscillator of base-pairs in the case of the strong coupling. The third point, the results of path 



integrals can be corrected by computing the higher terms in the expansion of the velocity equation (3.13) 
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